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DETECTION  OF  ARCTIC  WATER  SUPPLIES 
WITH  GEOPHYSICAL  TECHNIQUES 

S A Arcone,  A | Delanev  and  P V Sellmann 


INTRODUCTION 

(■roundwatri  distribution  patterns  in  arctic 
regions  are  similar  to.  but  more  restnc  ted  than 
those  in  geologically  analogous  permafrost  tree 
areas  (Hopkins  et  al  1‘tSS)  Where  permafrost  ex 
ists.  the  distribution  ot  groundwater  ran  be 
direr  fly  related  to  thaxs  /one  or  r urrenr  e whit  h 
depei  s on  the  physiography  setting  (slope,  ex 
posure,  vegetation,  soil  type,  proximity  to  maior 
water  bodies,  etc  ),  mean  annual  air  tempera 
true,  and  historical  fat  tors  However,  establish 
ing  the  distribution  ot  thaw  /ones  alone  does  not 
assure  a developable  source  ot  water  Ibe 
sediments  should  ideally  be  c oarse  grained  and 
extensive,  and  the  water  itxelt  should  be  tree  ot 
organic  and  mineral  pollution 

Situations  that  exhibit  a potential  for  water 
development  in  permafrost  areas  are  1)  un- 
frozen permeable  materials  above,  within  and  or 
below  the  permafrost.  2)  ponds,  lakes  and  risers 
that  do  not  tree/e  to  the  bottom  in  the  winter 
and  are  deep  enough  to  provide  a supply  ot 
water  without  high  concentrations  ot  impurities, 
and  t)  unfrozen  zones  below  river  and  lake  hot 
toms  where  a reservoir  of  water  may  be 
available  in  permeable  material  Most  of  these 
situations  represent  an  interruption  in  per 
mafrost  distribution  or  can  be  related  to  ice 
thickness  ot  water  bodies,  suggesting  that  their 
occurrence  should  decrease  as  one  moves  pro 
gressively  north 

there  are  many  geophysical  exploration 


methods  available  to  aid  in  locating  thaw  zones 


could  be  applicable  to  groundwater  surveys 
I hese  methods  identity  material  types  by  deter 
mining  physical  properties  such  as  acoustic 
velocities  (the  seismic  method)  or  electrical 
resistivity  (galvanic  electromagnetic  methods) 
I samples  ot  groundwater  investigations  in 
nonpermafrost  areas  using  these  methods  in 
elude  the  location  ot  buried  stream  channels 
determinination  ot  the  positron  of  the  water 
table  and  studies  ot  coastal  salt  water  intrusion 
Recent  surveys  ot  the  extensive  literature  and 
tec  hniques  utilized  tor  groundwater  detection  in 
both  permafrost  and  nonpermatrost  areas  have 
been  made  by  Zohdy  et  al  (1^741  and  Ogilvv 
(1^7) 

flee  tromagnetic  methods*  ot  geophysical  ex 
ploration  have  become  popular  since  the'  early 
to  mid  1%0\  because  ot  their  noncontact 
nature  (i  e lack  ot  earth-electrode  contact  and 
its  inherent  resistivity  problems)  and  their  well 
defint'd  depths  ot  penetration  These  methods 
utilize  propagating  elec  tromagnetic  energy  from 
the  entire  known  frequency  spec  trum  from  ultra 
low  frequency  signals  generated  by  lightning  to 
gamma  ray  radiation  produced  by  radioactive 
dec  ay  Methods  may  be  either  passive,  sue  h as  in 
the  detection  and  measurement  ot  telluric  cur 
rents,  self-potentials,  magnetic  fields  and  optic  al 
and  nuclear  radiation,  or  t hoy  may  be  active, 
such  as  in  the  transmitting  and  ree  owing  ot  radar 
signals  or  radiowaves  and  in  the  induction  and 
measurement  of  magnetic  fields  and  electric 
field  polarization  In  all  electromagnetic 


"Also  cellist  a i cm  alternating  > intent  mrihotts 


geophysical  methods,  the  physical  property 
measured  usually  depends  on  the  frequency 
used  General  information  on  most  methods  is 
given  in  Barnes  (1965),  Ogilvy  (1967),  Zohdy 
(1974).  and  Ward  (1967) 

The  propagation  of  electromagnetic  energy 
up  to  about  3X10'“  Hz  in  earth  materials 
depends  strongly  on  elec  trical  resistivity  (the  in- 
verse of  conductivity)  and  on  dielectric  proper- 
ties Propagation  also  depends  on  geometric 
considerations,  such  as  comparison  of 
wavelength  to  the  size  of  the  feature  being  in- 
vestigated Resistivity  and  dielei  trie  properties 
art'  strongly  influenced  by  the  presence  of  water, 
because  the  dissolved  ions  can  make  a conduc- 
tive, polarizable  electrolyte  Below  about  I 
MH/.  resistivity  is  the  predominant  property 
Resistivity  surveys  can  provide  information  on 
the  distribution  and  quality  of  subsurface  water, 
providing  that  resistivity  contrasts  exist  between 
the  water  or  water-bearing  strata  and  the  other 
surrounding  materials  Above  about  100  MH/. 
dielectric  properties  are  most  important 
Although  present  systems  (interferometric  and 
radar)  are  rarely  used  for  measuring  dielectric 
properties,  their  performance  depends  on  inter 
facial  dielectric  contrasts,  thus  enabling  con- 
tinuous profiling 

In  this  report  we  review  the  theory  and  ap- 
plication to  arctic  groundwater  problems  of 
several  active  techniques  used  in  the  radiowave 
band  between  15  kH/  and  10  CH/  The  studies 
presented  are  extracted  from  research  protects 
that  we  and  our  colleagues  at  CRREL  and  the 
Geological  Survey  of  Canada  have  performed 
over  the  past  several  years  The  discussions  ot 
data  interpretation  are  limited  to  brief  descrip- 
tions of  material  type  correlations  with  resistivi- 
ty readings  or  with  signal  strengths  More 
thorough  discussions  of  electromagnetic  field 
penetration,  interpretation  of  stratification,  and 
the  effects  of  topography,  lateral  inhomogeneity 
and  material  type  may  be  found  in  the 
references  cited  with  each  example 

The  body  of  this  report  is  divided  into  two 
parts  Part  I deals  with  low  frequency  (15  to 
400-kHz)  radiowave  and  magnetic  induction  re- 
sistivity methods  Part  II  deals  with  high  frequen- 
cy (0.1-  to  10-GHz)  radar  dielectric  methods 
Each  part  begins  with  a discussion  of  the  ap- 
propriate earth  electrical  properties,  followed 
by  a description  of  the  principles  and  equipment 
used  and  documented  studies  appropriate  to 
arctic  groundwater  investigation 


A review  ot  a CRRE  L study  ot  NASA  satellite 
imagery  over  arctic  lakes  is  presented  in  Appen- 
dix A Although  photographic  imagery  d - es  not 
directly  reveal  subsurface  information,  it  is  in 
eluded  because  of  its  application  to  arctic 
groundwater  problems 

PART  I.  LOW  FREQUENCY 
RESISTIVITY  METHODS 

Resistivity  of  earth  materials 

The  geophysical  systems  to  be  discussed  in 
this  section  measure  the  physical  property  of 
electrical  resistivity  at  frequencies  below  about 
1 MH/  We  therefore  begin  with  a discussion  of 
the  elec  trical  resistivity  of  earth  materials  with 
emphasis  on  materials  in  a frozen  siate  The  pro- 
perty of  dielectric  permittivity  has  a negligible 
effect  on  propagation  at  the  frequencies  con- 
sidered. while  the  property  of  magnetic 
permeability  is  essentially  the  same  in  most 
materials  throughout  nature,  the  maior  excep- 
tions being  a few  iron-bearing  minerals  More 
complete  discussions  of  electrical  properties  of 
earth  materials  may  be  found  in  the  texts  by- 
Keller  and  Erischknecht  (1966)  and  Parkhomenko 
(1967) 

The  electrical  resistivity  (or  its  inverse,  con- 
ductivity) of  a material  is  a bulk  property  ex- 
pressed in  ohm-meters  It  is  generally  inversely 
proportional  to  the  product  of  the  concentration 
of  c harge  carriers  and  their  mobility  Charge  car- 
riers are  termed  either  electronic  when  they  are 
electrons  (as  in  a metal),  or  electrolytic  when 
they  are  ions  (such  as  dissolved  salts)  The 
mobility  of  the  charge  carriers  is  determined  by 
many  factors,  such  as  temperature,  the  viscosity 
of  the  dissolving  fluid,  the  weight  of  the  ion,  or 
the  texture  of  the  material  through  which  the 
ionic  solution  must  pass 

Generally  natural  material  resistivity  is  deter- 
mined by  texture  (grain  size,  porosity  and 
permeability)  and  the  ionic  concentration  within 
the  adsorbed  water  on  grain  surfaces  and  within 
the  absorbed  water  in  pore  spaces  Electrical 
conduction  in  rocks  and  sediments  is  therefore 
mainly  electrolytic  and  takes  place  along  con- 
tinuous films  of  water  adsorbed  on  grain  bound- 
aries, as  well  as  through  water  that  may  fill  the 
available  pore  spaces  A major  mineralogical 
consideration  is  that  clays  can  have  a large 
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RESISTIVITY  (ohm-m«t«rf) 


Figure  1.  Resistivity  ranges  of  most  unfrozen 
earth  materials  (after  Culley  et  al.  1975). 


Figure  2.  Resistivities  of  several  saturated 
soil  types  and  one  rock  type  as  a function  of 
temperature  (compiled  by  Hoekstra  et  al. 
1974) 


amount  of  ions  adsorbed  on  their  crystal  lat- 
tices, producing  a surface  conduction  effect  that 
can  appreciably  lower  the  resistivity  of  a sedi- 
ment or  rock. 

In  Figure  1 the  expected  ranges  of  resistivity 
based  on  surveys  and  laboratory  measurements 
are  compiled  for  some  unfrozen  earth  materials. 
Among  the  unconsolidated  materials,  clay  and 
marl  through  loose  sands  generally  show  an  in- 
crease in  resistivity  with  increase  in  particle  size. 
This  may  be  expected  because  of  the  corre- 
sponding decrease  in  available  adsorbing  sur- 
face area  Alluvial  deposits  and  glacial  till  ex- 
hibit a wide  range  of  resistivities  due  to  large 
variations  in  textural  composition  The  con- 
solidated materials,  chalk  through  crystalline 
rocks,  generally  show  an  increase  in  resistivity 
with  decreasing  porosity  and  permeability  Clay 
content  may  also  be  significant  in  some  rock 
types  (Ward  1967),  and  the  large  amounts  of 
iron-  and  magnesium-rich  minerals  in  some 
basaltic  igneous  rocks  can  make  them  less 
resistive  than  granitic  igneous  types.  In  all  cases, 
the  effects  of  weathering,  jointing  or  fracture, 
which  increase  permeability  and  cause  the  for- 
mation of  clay  minerals,  significantly  lower  the 
resistivity  of  the  fresh  rock. 

In  frozen  soils,  temperature  and  ice  content 
greatly  influence  electrical  resistivity.  At 
temperatures  above  freezing,  only  small  in- 
creases in  resistivity  may  take  place  with 
decreasing  temperature.  Below  freezing, 
however,  temperature  changes  can  cause 
substantial  increases  in  resistivity,  as  is  il- 
lustrated in  Figure  2 for  several  saturated  sedi- 
ment types  and  one  rock  type.  Below  0°C,  the  in- 
crease in  resistivity  is  due  to  the  decrease  in  the 
amount  of  unfrozen  water  remaining  in  the 
ground  (Hoekstra  and  McNeill  1973,  Anderson 
and  Morgenstern  1973)  As  this  unfrozen  water  is 
in  the  adsorbed  state,  resistivity  generally 
decreases  with  decreasing  grain  size  (or  increas- 
ing area  of  surface  adsorption)  as  it  does  for  un- 
frozen sediments. 

Frozen  ground  may  also  have  ice  contents  in 
excess  of  saturation,  as  is  shown  by  the  occur- 
rence of  massive  ice  lenses  and  ice  wedges 
which  commonly  occur  in  fine-grained 
sediments  Hoekstra  et  al  (1974)  investigated  the 
relationship  of  resistivity  to  ice  volume  for  Fair- 
banks silt  and  found  an  approximate  linear 
dependency,  with  electrical  resistivity  varying 
from  about  2500  ohm-m  to  greater  than  100,000 
ohm-m  for  ice  volumes  between  60  and  100%. 
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figure  J.  Ground  electromagnetic  held  components  oi  a radio 
nave  transmitted  from  a vertically  polarized  antenna 


Theory  and  description  of  techniques 

Surface  impedance  hadiowavej  technique 

The  surface  impedance  technique  derives 
Kround  resistivity  values  from  a comparison  be- 
tween the  electromagnetic  field  components  of 
a propagating  ground  or  sky  wave  Between  15 
and  25  kHz  in  the  VIF  (very  low  frequency)  band, 
powerful  transmitters  operated  for  military  com- 
munications allow  radiation  to  be  monitored 
over  a range  of  several  thousand  kilometers  via 
the  sky  modes,  which  consist  of  waves  singly  or 
multiply  reflected  from  the  ionosphere  and  the 
earth  Station  NIK  (18b  kHz)  in  |im  Creek, 
Washington,  is  usually  monitored  throughout 
Alaska  Between  200  and  41 5 kHz,  numerous  low 
power  transmitters  operated  for  small  plane 
navigation  allow  radiation  to  be  monitored  only 
within  their  immediate  vicinity  This  radiation  is 
in  the  ground  wave  mode  wherein  the  radiation 
spreads  radially  over  the  ground  from  the  trans- 
mitter whose  range  is  dependent  on  ground 
resistivity 

The  electromagnetic  field  components  of  a 
ground  or  sky  wave  radiated  by  a vertically 
polarized  antenna  are  illustrated  in  Figure  3 £ 
refers  to  the  electric  field  components  and  H to 
the  magnetic  field  component  referenced  to 
local  x.  y,  z coordinates  The  surface  impedance 
Z.  is  then  defined  as* 

Z.*£./H,U  (1) 


*Ttm  definition  dutomdticdlly  applies  to  ihe  ground  wave, 
but  for  the  sky  wave  the  gracing  incidence  angle  upon  the 
earth  must  be  neat  0°.  which  is  usually  true  at  ranges  > UXX) 
km  at  Vlf 


Z.  is  a complex  quantity  having  both  an 
amplitude  |Z.|  and  a phase  i expressible  as  Z.  = 
|Z.|  e‘*  At  frequencies  below  about  1 MHz,  Z.  is 
approximated  for  a uniform  subsurface  model 
by 

Z.  = -v/cupop  e'4'*  (2) 

where  i = n/~-T 

g = resistivity  in  ohm-meters 
Mo  = free  space  permeability  (4n  X 10'T 
H m") 

to  •=  frequenc  y in  radians/second 

Equation  2 applies  when  dielectric  properties 
are  neglected  The  phase  value,  45°,  is  usually 
indicative  of  resistive  homogeneity,  at  least  to  a 
depth  6.  defined  in  the  literature  as  the  skin 
depth  (see  lackson  1%2  or  Stratton  1941)  and 
derived  from  the  formula 

d = V so  Mo  (3) 

At  this  depth,  the  refracted  field  strengths  have 
attenuated  to  e'1  of  their  surface  values 
Because  d depends  upon  frequency,  more  than 
one  frequency  can  be  used  to  resolve  resistivity 
changes  with  depth  In  Alaska,  adequate  field 
strengths  above  10  kHz  within  the  VLF.  LF  and 
MF  bands  are  rarely  simultaneously  available  at 
any  one  location 

When  the  earth  is  layered,  eq  2 must  be 
modified  because  the  phase  and  amplitude  of  Z, 
change,  with  the  magnitude  of  change  depend- 
ing upon  the  resistivity  and  thickness  of  each 
layer  Wait  (1962)  presents  formulas  for 
generating  Z,  values  above  any  number  of 
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layers  In  general,  a phase  greater  than  45° 
usually  indicates  that  resistivity  is  decreasing 
with  depth,  whereas  a phase  less  than  45°  usual- 
ly indicates  that  resistivity  is  increasing  with 
depth  According  to  this  simplified  model  of 
homogeneous  layers,  all  phase  angles  occur  be- 
tween 0°  and  90°  When  lateral  inhomogeneities 
exist,  the  theory  does  not  apply  and  only 
qualitative  interpretations  can  be  made 

The  instruments  used  for  measuring  resistivity 
in  the  VLF  and  LF  bands  are  similar  to  each  other 
and  commercially  available  E,  is  determined 
between  two  probes  inserted  in  the  ground  and 
compared  to  Ff,  measured  with  a ferrite-loaded 
coil  The  LF  device  must  be  pretuned  for  opera- 
tion between  200  and  400  kHz,  whereas  the  VLF 
device  operates  with  pretuned  circuits,  each 
pretuned  to  one  of  the  few  VLF  transmitters  in 
use  today  The  amplitude  of  Z,  is  calibrated  in 
ohm-meters  of  apparent  resistivity  e • (range  is 
0-30,000  ohm-m)  which  is  determined  from  the 
amplitude  inversion  of  eq  2 

0.  = |Z,|V«h.  (4) 

This  value  equals  the  true  resistivity  when  the 
earth  is  homogeneous  The  phase  of  Z,  is  read 
directly  in  degrees  (range  is  0°  to  90°)  A high  im- 
pedance voltmeter  is  used  to  measure  E,  so  that 
the  contact  resistance  of  the  probes  does  not  in- 
terfere with  the  measurement. 

The  accuracy  of  the  readings  depends  upon 
the  signal-to-noise  ratio.  This  ratio  is  mainly 
determined  by  transmitter  distance  and  power 
(for  LF  and  VLF),  by  ionospheric  factors  (for  VLF), 
and  by  the  value  of  ground  resistivity  itself  (f or 
LF  and  VLF)  because  as  resistivity  decreases  so 
does  the  amplitude  of  £„. 

Airborne  radiowave  technique 

The  surface  impedance  resistivity  method  has 
also  been  adapted  to  airborne  use  (Hoekstra  and 
McNeill  1973,  Hoekstra  et  al.  1974).  In  this 
method,  however,  E,  is  compared  to  E,  rather 
than  H,  (see  Fig.  3)  because  E,  cross-couples  dur- 
ing aircraft  roll  into  the  horizontal  antenna  that 
measures  £,.  Consequently,  only  the  component 
of  E,  that  is  out  of  phase  with  £,  can  be  used  £,  is 
then  used  as  the  phase  and  amplitude  reference 
against  which  £,  is  measured  so  as  to  avoid  mak- 
ing measurements  of  the  absolute  electric  field 
strength,  which  varies  with  range,  ionospheric 
factors,  etc.  (this  is  discussed  more  thoroughly  in 


Arcone  1977)  The  disadvantages,  then,  of  this 
system  are  that  only  a partial  measure  of  the  ap- 
parent resistivity  is  obtained  with  no  phase  infor- 
mation, and  that  tF*ere  is  a loss  of  ground  detail 
due  to  altitude  ettects  The  advantage  is  the 
greater  amount  of  area  that  can  be  covered 

Magnetic  induction  technique 

This  method  derives  ground  resistivity  values 
from  the  amount  of  magnetic  field  coupling  be- 
tween two  loop  antennas  located  at  or  slightly 
above  the  earth's  surface  One  loop,  the  trans- 
mitter, generates  a primary  a c.  magnetic  field 
that  couples  directly  to  the  receiver  loop 
through  free  space  and  also  induces  eddy  cur- 
rents within  the  earth  These  currents  then 
generate  a secondary  magnetic  field  that  also 
couples  with  the  receiver  and  is  used  in  com- 
parison with  the  primary  field  for  computing 
resistivity  values.  This  secondary  field  is  not  in 
phase  with  the  primary  field,  and  thus  it  has  both 
in-phase  and  quadrature  phase  components  with 
respect  to  the  primary  field  The  quadrature 
phase  component  is  used  for  computing  resistivi- 
ty values  because  it  has  a greater  range  of  linear 
response  with  resistivity.  In  addition  to  ground 
resistivity,  secondary  coupling  is  also  affected 
by  loop  separation  and  oriertation,  loop  height 
above  ground,  and  transmitter  frequency 

The  Geonics  EM-31  was  used  for  the  magnetic 
induction  surveys  This  instrument  is  calibrated 
in  millisiemens/meter  of  apparent  conductivity 
o.  (o.  = 1/p.)  At  the  instrument's  operating  fre- 
quency of  40  kHz,  the  apparent  conductivity 
with  respect  to  a first  layer  thickness  t for  a two- 
layer  model  varies  as  is  shown  in  Figure  4.  o.  is 
normalized  by  the  conductivity  of  the  first  layer 
and  t is  normalized  by  the  coil  separation  s = 
3.66  m As  shown  in  Figure  4,  when  a conductor 
lies  above  an  insulator  (o,/o,<1 ),  little  change  in 
o„  will  occur  for  values  of  t/s> 2 (or  f = 7.3  m for 
the  EM-31). 

Case  studies 

Geophysical  field  investigations  conducted  by 
CRREL  since  1971  have  been  performed  primari- 
ly in  permafrost  environments.  Some  investiga- 
tions have  been  conducted  in  temperate  areas  at 
the  request  of  other  agencies.  Many  aspects  of 
these  studies  are  applicable  to  arctic  ground- 
water  investigations  and  are  therefore  pre- 
sented The  studies  provide  information  on 
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Figure  4.  Two-layer  curves  for  a coplanar,  double  loop 
magnetic  induction  system  operating  at  a frequency  of  40 
kHz.  For  the  Ceonics  EM-31,  s = 3. 66  m.  ( Courtesy  of 
Ceonics  Ltd.) 


water  below  lake  and  river  ice,  thaw  zones 
below  the  beds  of  streams  and  lakes,  unfrozen 
zones  in  permafrost  areas,  and  distribution  of 
materials  that  have  high  permeability  and  may 
provide  good  groundwater  sources. 


Location  of  thaw  zones  beneath  lakes  and  rivers 

Smith  Lake,  Fairbanks,  Alaska.  Smith  Lake  is 
about  10  km  west  of  Fairbanks  in  the  discon- 
tinuous permafrost  zone.  The  permafrost,  a 
frozen  silt,  surrounds  the  lake  (P4w4  1958).  At 
the  time  of  the  surveys,  lake  ice  thickness  varied 
between  66  and  76  cm  and  water  depth 
increased  gradually  from  the  margins  to  the 
center  where  the  depth  was  approximately  2 m 
including  the  ice. 

The  lake  was  investigated  in  April  1974  using 
the  surface  impedance  method  at  18.6  kHz  and 
in  April  1977  using  the  magnetic  induction 
(EM-31)  method.  Over  the  central  portion  of  the 
lake,  two  VLF  readings  gave  110  and  120  ohm-m 
at  45°  (very  few  readings  were  successful  due  to 
the  weak  reception  in  the  area)  while  the 
magnetic  induction  readings  made  at  3-m  inter- 
vals (Fig.  5)  ranged  between  130‘and  180  ohm-m. 
Over  the  permafrost,  one  VLF  measurement  was 


read  at  between  500  and  1000  ohm-m  at  about 
45°,  while  the  EM-31  readings  ranged  from  400 
to  700  ohm-m. 

The  VLF  data  imply  deep  thaw  under  the  lake 
and  probably  no  permafrost.  At  120  ohm-m,  the 
skin  depth  at  18.6  kHz  is  40  m which  is  close  to 
the  maximum  permafrost  thickness  in  other 
parts  of  this  area  The  magnetic  induction 
resistivities  are  greater  here  because  of  the 
greater  sensitivity  to  the  ice  cover  and  because 
the  instrument  was  held  at  a 1-m  height.  This  is 
verified  in  Figure  6 which  shows  that  a model 
120-ohm-m  real  resistivity  for  the  water  and 
thawed  bottom  sediments  rises  to  a 165-ohm-m 
apparent  resistivity  with  an  ice  cover,  which  is 
typically  about  bOOO  ohm-m.  Increases  in  ice 
resistivity  above  about  1000  ohm-m  marginally 
affect  the  readings. 

Although  Smith  Lake  has  over  1 m of  unfrozen 
water  remaining  near  its  center  during  late 
winter,  the  resistivity  data  do  not  indicate  its 
presence  in  spite  of  giving  information  on  depths 
of  thaw  This  is  due  to  the  lack  of  resistivity  con- 
trast between  the  water  and  sediments  or  even 
within  the  sediments  themselves. 
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fillin'  S Magnetic  induction  apparent 
resistivity  traverse  over  Smith  Lake,  near  lair 
hanks.  Alaska 


i ?,  lohm-m  \ W! 


figure  t>.  Magnetic  induction  apparent  resistivity 
as  a function  of  ice  resistivity  for  a lake  model  in 
which  wafer  and  bottom  sediments  are  of  the 
same  resistivity  Resistivity  variations  above  IOOO 
ohm  m in  the  ice  cover  have  marginal  effect  upon 
the  apfrarent  resistivity 
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Thaw  Takes.  near  Pump  Station  S.  Aheska 
Pipeline  Route.  Alaska.  1 he  results  of  this  study 
(Sellm.rrm  et  al  1474)  are  shown  in  ft  (jure  7 
Resistivity  profiles  (f  ig  7h  and  c)  using  the  stir 
face  impedance  ter  hmque  at  VI  f were  made 
across  two  small  lakes  near  Prospect  C amp  on 
the  Trans  Alaska  Pipeline  route  near  pump  sta 
lion  S (fig  7a)  These  are  prohahly  thaw  lakes 
surrounded  by  permafrost  Apparent  resistivity 
readings  over  the  surrounding  tundra  average 
about  TOO  ohm  m,  but  these  values  drop  to  as 
low  as  t S ohm  m along  sei  tion  a-a'  and  as  low  as 
SO  ohm  in  along  section  I B over  the  lakes  No 
investigation  was  made  of  the  amount  ot  water 
under  the  ice  or  of  the  depth  of  thaw  beneath 
the  bed  However,  the  very  low  resistivities 
measured  along  a-a'  are  characteristic  of  clay 
and  silt,  implying  that  some  thaw  exists 
Therefore,  it  water  is  present  under  the  ice,  it 
probably  occurs  in  thawed,  finegrained 
sediments  trom  which  water  would  be  difficult 
to  extract 

Lake  bed.  Culkana.  Alaska  This  study  (Arc  cine 
et  al  1‘>78)  was  conducted  at  a small  lake 
located  near  the  town  of  Culkana  (fig  8)  The 
soil  types  of  the  area  are  generally  described  as 
glacial  clays,  day  tills  and  lacustrine  deposits 
that  arc>  fairly  homogeneous  in  nature  exc  ept  tor 
their  ice  content  lho  site  is  situated  along  one 
of  the  earlier  proposed  Trans  Alaska  Pipeline 
alignments  An  interpretation  of  2 S borings 
(R&M  f ngineering  WO)  from  this  alignment  is 
presented  in  f igure  8b 

figure  8a  shows  the  apparent  resistivity  pro- 
files measured  by  magnetic  m<  ac  tion  in  luly  of 
f‘)7b  when  the  active  layer  was  thawed  and 
again  in  April  of  1 **77  when  the  active  layer  was 
frozen  The  qualitative  nature  of  both  resistivity 
profiles  remains  the  same,  but  the  greater 
penetration  afforded  in  April  (when  the  ac  five 
layer  was  frozen)  allowed  more  sensitivity  to 
deeper  soil  variations,  especially  for  the  higher 
icc  , ontent  regions  I he  ettec  t of  the  large  thaw 
zone  at  the  center  is  to  lower  the  apparent 
resistivity  to  less  than  SO  ohm  ni  during  both 
observation  pen**!'  I his  low  value  of  resistivity 
is  indie  alive  of  the  tine  grained  sediments  of  the 
thaw  zone  and  probably  indicates  a poor  water 
supply  It  is  similar  to  the  thaw  zone  disc  ussed 
previously 

Thaw  lake.  Alaskan  Arctic  Coastal  Plain  A 
comparative  survey  was  conducted  at  an  ar 
tific tally  deepened  lake  near  Prudhoe  Hay. 
Alaska,  situated  next  to  the  Sagavamrktok  River 


Resistivity  was  measured  using  the*  l 7 m coil 
spat  mg  magnetic  indue  lion  instrument  and  the 
surf  ac  e impedanc  e method  at  2 14  kHz  * igure  '*a 
shows  the  data  profiles  and  the  t ross  set  lion  ot 
lhi>  lake'  beneath  the  profile's  Also  shown  in 
I igure  Ma  is  the  surfat  e imped  am  e phase  profile 

I igure  ‘lb  shows  a vertic  ally  layered  resistivity 
model  ot  the  ARl  O lake  and  its  subsurface 
sediments  with  the  actual  values  ot  n e and 
wafer  thicknesses  Theoretical  surface  im 
pedam i‘  and  magnetic  induction  vain  s,  which 
are  similar  to  the  actual  data  taken  ovet  the 
lake,  arc*  given  above  the  model  A three  layer 
c ase  is  given  here  but  other  models  may  also 
produce  the  same  data  I he  upper  layer  is 
relatively  pure  freshwater  ice  with  resistivities 
usually  nun  h greater  than  the  value  given 
However,  lor  this  t ase  ot  an  insulator  over  a t on 
due  tor,  both  tec  hniques  aie  insensitive  to  the 
resistivity  value  ot  the  insulator  when  it  is  10,000 
ohm-m  or  greater  hut  are  sensitive  to  the'  ice 
depth  1 he  second  layer  is  ot  fresh  water  but  at 
an  unusually  low  resistivity  value  ot  40  ohm  m. 
suggesting  that  the  water  has  a very  high  ion 
content  I he  third  layer  is  thawed  sediment 
which  may  be  as  deep  as  20  to  TO  m I his 
material  is  composed  primarily  of  sands  and 
gravels  tor  which  1000  ohm  m is  not  unreason 
able  In  this  case,  the  value  of  1000  ohm  m was 
required  to  bring  the  model  phase  data  well  into 
the  observed  range  Therefore,  this  thawed 
material  is  most  likely  highly  permeable  and 
capable  ot  holding  a great  deal  ot  water  No  in 
termediate  layer  ot  finegrained  lacustrine 
sediments  is  shown  at  the  lake  hed  in  the  model 
As  this  is  ail  artificially  deepened  lake,  such  a 
layer  would  be  ot  only  marginal  thickness 

lhe  large  degree  of  variability  in  the  II  data 
relative  to  the  magnetic  indue  tion  data  deserves 
some  comment  Near  the  shore  the  apparent  re- 
sistivity at  2T4  kHz  shows  a very  interesting  in- 
terference pattern  whir  h is  usually  indicative  ot 
a subsurfac  e lateral  inhoniogeneity  in  resistivity 
Over  the  central  portion  ot  the  lake,  however, 
the  source  of  the  variability  is  not  certain  One 
possibility  is  an  uneven  distribution  in  thickness 
and/or  resistivity  of  the  bottom  sediments  Note 
that  the  lowest  values  of  resistivity  are 
associated  with  the  highest  values  ol  phase, 
which  is  thc'oretic ally  correct  in  view  ot  the 
resistivity  stratification  (see.  tor  example. 
Arconeetal  WB) 

lhassed  rivet  channels  Sagavanirktok  and 
kuparuk  Risers  Small  water  bodies  can  uwially 


a.  Location  map  showing  survey  lines  near  pump  sta-  b.  Apparent  resistivity  profile  at  VLF  along  line  E-B 
lion  5. 


Dittonc*,  rnliii 


c.  Idealized  section  and  associated  apparent 
resistivities  across  "Extra  Lake"  along  line  a-a'. 


Figure  7.  Resistivity  profiles  using  the  surface  im- 
pedance technique  of  two  thaw  lakes  near  pump  sta- 
tion 5,  Alyeska  Pipeline  route,  Alaska. 


a Resistivity  profiles  b Geologic  cross  section 


figure  H Magnetic  induction  apparent  resistivity  profiles  lal  taken  user  a small  thass  lake  near 
Culkana.  Alaska  Corresponding  geologic  cross  section  lb)  alter  R\M  l ngmeering  and  Ceological 
Consultants  ( 1970) 


b Layered  geoelectric  model  of  the  ARCO 
a.  Surface  impedance  apparent  resistivity  Lake  showing  layer  depths,  resistivities  and 

and  phase  and  magnetic  induction  aoparent  dielectric  permittivities  The  permittivity 

resistivity  profiles  across  the  ARCO  Lake.  values  (k)  are  of  little  significance  at  the  tre- 

(Approxima te  lake  cross  section  along  quences  used.  The  values  given  above  the 

traverse  from  data  courtesy  of  Austin  model  are  theoretical  and  are  to  be  compared 

Kovacs,  CRREL).  with  the  actual  data  of  9a. 

figure  9.  Apparent  resistivity  profiles  and  layered  geoelectric  model  of  ARCO  l ake  at  Prudhoe  Hay. 

Alaska. 
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figure  10.  VLF  apparent  resistivity  profile  across  the  main  channel  of  the 
Sagavanirktok  River  near  Franklin  Bluffs.  Alaska 


be  found  in  the  winter  only  beneath  major  arctic 
streams  because  most  streams  freeze  to  their 
bed  Thus,  seasonal  discharge  of  arctic  rivers 
and  streams  usually  makes  for  an  unreliable 
winter  water  supply  In  some  locations,  however, 
the  alluvium  beneath  the  bed  can  be  thawed  to 
some  depth,  forming  small  reservoirs  and 
groundwater  aquifers.  Geophysical  techniques 
may  then  be  used  to  detect  such  thaw  zones 
since  they  are  more  conductive  than  the  sur- 
rounding material 

Figure  10  (after  Delaney  et  al.  1975)  shows  a 
VLF  (186-kHz)  apparent  resistivity  profile  per- 
formed in  April  1974  across  the  Sagavanirktok 


Figure  1 1.  Magnetic  induction  profile  across 
the  Kuparuk  River  (British  Petroleum  Co. 
water  supply],  near  Prudhoe  Bay,  Alaska 
(after  Northern  Engineering  Services  Com- 
pany 1977). 


River  near  Franklin  Bluffs  in  northern  Alaska 
The  deeper  and  more  active  channel  positions 
correspond  to  the  lower  resistivity  values  The 
lowest  values,  for  the  channel  on  the  left,  may 
suggest  deeper  thaw. 

A similar  study  (Fig.  11)  using  the  magnetic  in- 
duction technique  was  conducted  by  Northern 
Engineering  Services  Company  (1977)  and  sup- 
ported by  a drilling  program  The  Kuparuk  River 
near  Prudhoe  Bay,  Alaska,  was  profiled 
Resistivity  minima  were  recorded  at  stations  50 
and  60  At  station  50,  a drill  hole  revealed  1.8  m 
of  ice  over  2.6  m of  water  At  the  end  stations  the 
ice  was  frozen  to  the  bed 
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figure  \2.  Apparent  resistivity  contours  at  17  8 kHz  over  a 12-acre  held 


in  Plainfield.  New  Hampshire. 


Location  of  permeable  materials 
in  unfrozen  ground 

A VIF  surface  impedance  survey  was  con- 
ducted in  Plainfield,  New  Hampshire  (Delaney  et 
al  1974).  to  locate  permeable  material  that 
might  be  recharged  from  a nearby  stream  The 
area  surveyed  was  approximately  12  acres 
covered  by  low  relief  glacial  drift  Resistivity 
contours  based  on  a 25-m  grid  data  pattern  are 
shown  in  Figure  12  Most  of  the  area  was 
relatively  uniform  in  resistivity  with  values  in  the 
200-400  ohm-m  range  On  the  eastern  margin  of 
the  study  area  adjacent  to  the  stream,  an  area  of 
higher  resistivity  between  500  and  900  ohm-m 
occurred 

Based  on  the  range  of  resistivities  for  various 


earth  materials  shown  in  Figure  1,  values  in  ex- 
cess of  500  ohm-m  in  this  type  of  setting  can  be 
associated  with  sand  and  gravel,  materials  of 
good  permeability  Two  test  pits  dug  with  a 
backhoe  at  stations  N and  S revealed  clean 
gravel  with  a high  groundwater  table  A test  hole 
drilled  at  station  D revealed  less  permeable 
material 

Delineating  permafrost  boundartes 

Fox,  Alaska  Fairbanks  and  its  outlying  areas 
lie  within  the  discontinuous  permafrost  rone  of 
Alaska  This  area  is  generally  mantled  with  silt 
that  overlies  a Precambrian  schist  known  as 
Birch  Creek  schist  Although  the  fine-grained 
sediments  in  this  region  do  not  have  a potential 
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Figure  IJ.  Apparent  resistivity  profile  with  magnetic  induction  along  a section 
of  the  Fairbanks-Fox  Road  centerline,  and  the  geological  cross  section  derived 
from  drill  logs  I numbered  in  the  figure I provided  by  the  Alaska  Highway  Depart- 
ment. iQf  — unfrozen  silt,  Qsu  — frozen  silt.  PC  be  — Birch  Creek  schist 
bedrock.) 


for  ground  water  development,  the  region  pro- 
vides many  good  examples  of  how  thaw  zones  in 
perennially  frozen  sediments  can  be  detected. 

A 610-m  linear  section  (Fig  13)  of  the  pro- 
posed Fairbanks-Fox  Road  near  Engineer  Creek 
was  surveyed  (Delaney  et  al  1977)  for  massive 
ice  using  the  magnetic  induction  technique 
Measurements  were  made  at  3-m  intervals  dur- 
ing early  April  when  the  surface  was  snow- 
covered  and  the  active  layer  was  frozen  to  its 
maximum  depth.  This  section  is  on  a north- 
facing hillside  and  encounters  a transition  from 
the  unfrozen  (Qf)  to  the  perennially  frozen  (Qsu) 
silt  units  The  depth  to  bedrock  increases  fairly 
uniformly  from  6 7 to  16  0 m and  several  zones 
of  massive  ice  exist  within  the  Qsu  unit  (borings 
by  the  Alaska  Highway  Department  1973, 1974). 
The  geologic  cross  section  was  derived  from  the 
borings  and  the  surficial  geologic  map  of  P^w^ 
(1958). 


From  0 to  1 20  m in  the  Qf  unit  all  readings  fall 
between  200  and  400  ohm-m,  and  between  120 
and  280  m in  the  Qsu  unit  the  readings  range 
from  200  to  1000  ohm-m  Previous  measure- 
ments made  by  Hoekstra  et  al  (1974)  in  the  Fair- 
banks area  using  the  galvanic  resistivity  method 
showed  similar  values  for  Qf  and  Qsu  materials. 
Beyond  280  m in  the  Qsu  unit,  the  readings  in- 
crease to  between  2000  and  10,000  ohm-m  This 
last  subsection  was  found  to  contain  massive 
ice,  the  distribution  of  which  is  inferred  from  the 
drill  logs.  The  study  therefore  demonstrates  not 
only  the  use  of  this  technique  for  detecting  per- 
mafrost margins  but  also  for  detecting  material 
variations  within  permafrost 

Tuktoyaktuk,  Canadian  NWT.  In  April  1973,  an 
airborne  VLF  survey  was  performed  near 
Tuktoyaktuk,  NWT,  Canada,  for  the  Geological 
Survey  of  Canada  (Scott  and  Hunter  1976)  This 
region  is  predominantly  flat  and  underlain  by 
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a.  Thawed  lake 


b.  Permafrost 


Figure  14.  VLF  apparent  resistivity  contours  in  the  Tuktoyaktuk  region  of  the  Canadian  Northwest  Ter- 
ritories above  a thawed  lake  lal  and  permafrost  fb).  Map  courtesy  of  Geological  Survey  of  Canada  Con- 
tour values  are  in  ohm-m , shaded  prot ions  are  areas  of  deep  thaw.  The  heavy  bar  in  b is  where  ground  VLF 
reconnaissance  was  made. 


ice-rich  clay  Jilt  lhaw  lakes  occurring 
throughout  the  region  have  had  a significant  in- 
fluence on  development  ot  the  present  land- 
scape Surface  impedance  investigations  by 
Scott  and  Hunter  showed  areas  of  permafrost  to 
have  VLf  apparent  resistivities  between  1000 
and  20,000  ohm  m.  but  no  VI  f measurements 
were  made  over  thaw  z ones 

A simplified  map  ot  contoured  airborne  data 
resistivity  superimposed  upon  shaded  areas  o< 
copied  by  thaw  lakes  is  presented  in  figure  14 
I he  steep  gradient  in  the  contours  ot  I igure  14a 
parallels  the  lake  margin  and  thus  enc  loses  the 
thaw  / one  caused  bv  the*  lake  figure  14b  shows 
an  area  iust  south  ot  the  area  ot  14a  that  is 
predominantly  tro/en  and  contains  much  ice 
near  the  surface  In  the'  lower  right  band  corner 
ot  figure*  14b  there  occurs  a thaw  /one  \sith  a 
relatively  low,  HXX)  ohm  m contour  lying  chrec  t 
ly  over  it  Ihus.  it  is  seen  that  thaw  /ones 
assoc  lated  with  contemporary  lakes  and  drained 
lake  basins  may  be  delineated  on  a large*  sc  ale* 
with  airborne  adaptations  ot  the  surface  im 
pedanc  e resistiv  ity  method 

PART  II.  HIGH  FREQUENCY 
DIELECTRIC  METHODS 

Dielectric  properties  of  earth  materials 

The  geophysical  exploration  systems  dis- 
cussed in  this  section  are  radar  systems  which 
are  used  for  continuously  profiling  or  mapping 
surface  and  subsurface  variations  in  electrical 
properties  of  earth  materials  As  with  surface  or 
dinance  radar,  the  backscatter  of  beamed 
energy  at  frequencies  usually  between  KXT  and 
lO.tXXl  MHz  is  monitored  Subsurface  reflections 
take  place  where  strong  contrasts  in  bulk  dielec- 
tric permittivity  or  in  bulk  conductivity  exist  At 
radar  frequencies,  the  dielectric  permittivity  is 
usually  more  important 

The  dielectric  permittivity  t is  a measure  of 
the  electrical  polarizability  of  a material  and 
determines  the  velocity  of  propagation  v in  that 
material  by  the  formula 

V = 

where  p is  the  magnetic  permeability  of  the 
medium  which  is  the  same  as  p«  of  eq  2 for  most 
materials  Normally,  c is  expressed  as 

* = **«.  (til 


where  x*  is  the  complex,  relative  dielectric  con 
stant  that  is  dimensionless,  and  c„  is  the  dielec 
tru  permittivity  of  tree  space  = 8 ttS  X 10  " 
f/m  These  units  may  also  be  expressed  as 
(couloumb  motors, meter1  per  volts  meter)  or 
dipole  moment  density  per  unit  of  elec  tric  held 
strength 

x*  is  a c omplex  quantity  bee  ause  the  polariza- 
tion of  a material  may  not  always  be  in  phase 
with  the  impressed  electric  field  In  such  situa- 
tions x*  is  frequency-dependent  and  usually  ex 
pressed  as 

x*  = x’  — rx*  (7) 

where  x.’  is  the  in-phase  or  real  part  of  the  dielec 
trie  constant  and  x*  is  the  quadrature  phase  or 
imaginary  part  A finite  value  ot  x‘  indicates  that 
electromagnetic  energy  will  be  lost  during  pro- 
pagation lhe  conductivity  o also  causes  at- 
tenuation and  the  two  losses  cannot  be 
separated  by  measurement  It  is  then  customary 
to  define  an  effective  quadrature  phase  term  as 
o.„  so  that 

x*  = x"  — ia,„  (fi) 

where 

o.„  = x"  + o/cet„ 

There  are  three*  ways  an  individual  material 
may  be  polarized  first,  it  may  have  the  centers 
ot  the*  negatively  charged  electron  clouds  ot  its 
individual  atoms  or  molec  ules  slightly  displaced 
relative  to  the  positively  c harged  nuc  lens  Sec 
ond.  a net  displacement  between  negative  and 
positive  ions  may  occur  Ihird  the  molecules 
may  be  intrinsically  polarized  but  oriented  ran 
domly  as  with  water  Therefore  an  impressed 
elec  tric  field  polarizes  the  material  by  giving  the 
molecules  a preferred  orientation 

In  composite  materials,  interfacial  polari/a 
tion  may  occur  when  the  charge  flow  is  impeded 
at  an  interface  Such  phenomena  are  usually 
confined  fo  lower  frequencies  than  those  con- 
sidered in  this  section  of  the  report 
. If  ts  custom  ay  to  present  the*  dispersion  or  fre- 
quency dependence  ot  the  complex  values  ot  x 
when  discussing  dielectric  properties  Because 
water  is  highly  dispersive  above  1(XT  MHz.  all 
water-bearing  earth  materials  show  some  degree 
of  dispersion  figure  IS  shows  the  dispersion  of 
x'  and  x*  for  water  The  dispersion  for  unfrozen 
water  results  from  the  molecular  orientation  of 


IS 


Figure  1b.  The  complex  dielectric  constant  of  Goodrich 
clay  at  a water  content  of  10%  (g  H,0/g  soil) as  a function 
of  frequency  at  two  temperatures  (°C). 


the  polar  water  molecule,  whereas  the  disper- 
sion for  ice  results  from  proton  migration.  Both 
types  of  dispersion  are  described  as  relaxation 
polarizations  which  means  that  the  polarization 
of  the  substance  can  no  longer  keep  in  phase 
with  the  applied  electric  field  near  and  above 
some  characteristic  frequency. 

As  shown  in  Figure  15,  the  dispersion  of  water 
mainly  takes  place  above  about  200  MHz  and 
reaches  a maximum  loss  (i.e.  maximum  value  of 
k’)  at  about  20  CHz.  Below  about  200  MHz 
water  is  essentially  lossless  with  a very  high 
value  of  «'  at  about  78  The  dispersion  of  ice 


takes  place  at  very  low  frequencies  Above 
about  01  MHz.  ice  behaves  like  a lossless  dielec- 
tric with  k'  equal  to  about  3. 

Figures  16, 17  and  18  show  the  dispersions  of  a 
clay  and  a silty  clay  at  different  temperatures 
and  water  contents  (after  Hoekstra  and  Delaney 
1974).  When  these  curves  are  compared  with 
those  of  Figure  15,  it  is  evident  that  the  range  of 
frequencies  over  which  the  dispersion  occurs  is 
lowered  from  that  of  unfrozen  water  by  factors 
between  about  5 and  20.  and  the  low  frequency 
asymptotic  values  of  k‘  and  it'  now  are  between 
that  of  water  and  of  pure  sediment  («’  is  about  3 
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Figure  17.  The  complex  dielectric  constant  of  Goodrich  clay 
at  20°C  as  a function  of  frequency  at  two  wafer  contents 
fg  HtO/g  soil) 
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Figure  18.  The  complex  dielectric  constant  of  Suffield 
silty  clay  at  a water  content  of  10%  (g  H,0/g  soil).  The 
solid  lines  are  curves  generated  by  fitting  two  Debye  type 
dispersions  to  the  experimental  points. 
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/ inure  19  Ground-based  impulse  radar.  I he 
dark  shaded  area  is  the  transmutes!  beam  and 
the  hnhter  shading  is  the  reflected  beam  I he 
line  of  impulse  waveforms  represents  the  pot 
ti on  of  the  beam  that  is  received  lime  of 


to  4 tor  dry  silt)  Note  that  at  -10°C  in  t igures  1b 
and  18b  there  is  appreciable  dispersion  which  is 
due  to  the  unfrozen  adsorbed  water  These 
curves  suggest  that  serious  attenuation  will  take 
place  tor  radar  signals  propagating  through 
frozen  sediments  between  0 1 and  10  GHz 

Dry  materials  are  generally  nondispersive,  but 
considerable  information  is  lacking  on  their 
behavior  at  microwave  frequencies  Often,  the 
constituent  minerals  show  some  degree  of 
anisotropy  At  frequencies  between  TOO  Hz  and 
10  MHz,  Parkhomenko  (19b7)  reports  that  values 
of  *'  most  generally  fall  between  5 and  10  for  dry 
sedimentary  rocks  with  the  exception  of  some 
shales  which  are  rocks  high  in  clay  mineral  con 
tent  The  addition  of  water  to  the  more  porous 
sedimentary  rocks  may  then  cause  «'  to  increase 
to  as  high  as  40  (Parkhomenko  19t>7) 

Theory  and  description  of  equipment 

Profiling  radar 

Profiling  radars  have  been  described  in  the  ar 
tides  by  Moftitt  and  Puskar  (197b).  Davis  et  al 
(197b)  and  Annan  and  Davis  (197b)  The  principal 
features  and  operation  of  a subsurface  profiling 
radar  system  are  shown  in  Figure  19  A transmit- 
ting and  receiving  antenna  can  be  mounted  on  a 
sled  and  continuously  towed  over  the  ground  or 
ice  surface,  or  it  can  be  attached  to  a slow  flying 
aircraft  The  transmitter  radiates  a continuous 
tram  of  pulses  that  are  monitored  by  the  receiver 
after  reflecting  from  subsurface  interfaces  The 
time  and  amplitude  of  return  are  graphically 
displayed  to  produce  an  apparent  profile  of  the 
subsurface  interfaces  If  the  dielectric  permit- 
tivity of  the  media  above  the  interface  is  known, 
delay  time  can  be  converted  into  distance 

In  practice,  the  impulse  radiated  is  about  one 
tuli  sinusoid  lasting  about  10  to  20  ns  The  fre- 
quency spectrum  is  centered  at  about  100  MHz, 
which,  as  was  shown  in  Figures  lb-18,  is  about 
optimum  for  minimal  ground  attenuation  The 
bandwidth  is  about  1 50  MHz  Recently,  sinusoid- 
al pulses  centered  at  about  b25  MHz  have  been 
used  over  sea  ice  (Kovacs  and  Morey  1978) 
These  pulses  have  bandwidths  of  about  500 
MHz,  which  allow  serious  pulse  distortion  when 
propagation  is  in  a dispersive  material  There- 
fore, this  frequency  is  best  applied  in  areas  con- 
taining ice  or  material  that  is  frozen  or  of  low 
moisture  content 

Ideally,  the  ground  surface  as  well  as  the  sub- 
surface interfaces  should  be  smooth  to  avoid 


noise  caused  by  mechanical  vibrations  and  un- 
wanted signal  backscatter  (termed  clutter)  Con- 
sequently, best  results  are  often  obtained  over 
smooth  ice  covers 

Imaging  radar 

Surveillance  radars  are  those  that  are  able  to 
construct  two-dimensional  images  from  the  time 
delay  and  amplitude  of  the  signal  return  In 
geological  reconnaissance,  a type  of  surveil- 
lance radar  used  is  known  as  STAR  or  side- 
looking airborne  radar,  illustrated  in  Figure  20 
In  a SIAR  system,  an  antenna  fixed  beneath  an 
aircraft  radiates  a beam  which  is  narrow  in  the 
direction  of  the  line  of  flight  and  broad  in  the 
transverse  direction  The  time  of  return  is  simply 
related  to  the  ground  distance  and  the 
amplitude  is  mainly  determined  by  the  orienta- 
tion of  the  local  ground  profile  with  respect  to 
the  beam  The  width  of  illuminated  terrain  is 
reasonably  fixed  by  the  beam  shape,  and  the 
scanning  swaths  are  parallel  w ith  the  flight  line 
These  systems  normally  operate  in  the 
microwave  band  between  1 and  10  GHz  The 
radiated  signals  are  a series  of  pulses  containing 
several  hundred  or  thousand  cycles  of  the 
microwave  frequency  Consequently,  beam- 
width  and  bandwidth  are  well-defined  Since 
earth  materials,  even  when  frozen,  can  exhibit 
appreciable  loss  at  these  frequencies,  most  sub- 
surface penetration  will  occur  only  through 
relatively  lossless  dielectrics  like  freshwater  ice. 
dry  sediment,  or  snow 

Case  studies 

Radar  profile  of  a river  channel 

The  Geological  Survey  of  Canada  (Annan  and 
Davis  1977)  has  used  the  impulse  radar  tech- 
nique to  profile  ice  depth  and  bathymetry 
beneath  an  ice  road  (or  bridge)  across  the  Yukon 
River  at  Dawson  City,  Y T , Canada  Two  dif- 
ferent profiles  were  made,  one  with  an  impulse 
of  2-ns  duration  and  one  with  a 10-ns  duration 
Figure  21  shows  the  2-ns  and  10-ns  profiles  and  a 
geological  interpretation  The  closely  spaced 
dark  bands  correspond  to  the  successive  max- 
ima of  a single  impulse  return  The  thin  white 
lines  are  the  points  of  zero  amplitude  within  the 
sinusoidal  pulse  The  various  interfaces  with 
which  the  returns  correspond  are  labeled  in  the 
figure 

In  comparing  the  2-ns  to  the  10-ns  profile,  it  is 
apparent  that  the  2-ns  system  does  not  have  the 


TRAVEL  TIME  (r.$) 


Figure  20.  Side-looking  airborne  radar.  Time 
of  return  is  translated  into  distance  and 
amplitude  of  return  into  image  contrasts. 
Major  radiation  lobes  are  shaded.  Return 
energy  is  backscattered  rather  than 
specularly  reflected. 


a.  Radar  section  of  the  ice  bridge  across  the  Yukon  River  obtained  with  the  2-ns  impulse  antennas. 


Figure  1.  Impulse  radar  profile  at  2-ns  and  10-ns  pulse  duration  across  the  Yukon  River  at  Dawson 
City,  Y.  T.  (courtesy  of  the  Geological  Survey  of  Canada). 


20 


0 

M 


lit  M 


SurfOC* 

Reflection 

Bottom 

Reflection 


Figure  22  Radar  profile  obtained  from  a helicopter  over  Big  Lake  on  the  Alaskan 
North  Slope  ( after  Kovacs  IB78)  The  bright  spots."  dark  areas  on  the  record,  in- 
dicate potential  sites  where  wafer  still  exists  under  the  ice. 


water  penetration  capability  of  the  10-ns  system 
The  2-ns  system  has  a center  frequency  of  about 
500  MHz  for  which  the  loss  factor  k"  (Fig  15)  ot 
water  has  a value  of  approximately  1 At  10  ns. 
the  center  frequency  is  about  100  MHz  for  which 
*'  of  water  is  negligible  Therefore,  a 10-ns  im- 
pulse system  is  more  ideally  suited  for 
freshwater  probing 

Impulse  radar  profile  of  a freshwater  lake 

CRREl  (Kovacs  1978)  has  made  radar  profiles 
from  a helicopter  of  Big  Lake  in  the  National 
Petroleum  Reserve-Alaska,  on  the  Arctic  Coastal 
Plain  The  lake  was  almost  completely  frozen  at 
the  time  of  the  survey  The  radar  profile  of 
Figure  22  shows  two  small  areas  of  strong  returns 
at  a two-way  time  return  of  24  ns  This  time 
return  corresponds  to  a depth  of  2 m,  consider- 
ing a value  of  3 2 for  k'  in  the  ice.  which  is  close 
to  the  actual  lake  depth  A drill  hole  at  the 
return  marked  "b"  on  the  figure  verified  the 
presence  of  4 cm  of  water  between  the  ice  and 
the  lake  bed 

SLAR  imagery  of  arctic  lakes 
During  the  spring  of  1974,  SLAR  imagery  was 
obtained  for  part  of  the  Arctic  Coastal  Plain 
eastward  along  the  coast  from  Pt  Barrow. 
Alaska  (Weeks  et  al  1977).  This  area  is  generally 
flat  and  covered  with  numerous  shallow  lakes  of 
depths  ranging  from  a few  to  tens  of  meters  At 


the  time  of  the  survey,  this  area  was  still  snow 
covered  and  thaw  lakes  had  their  maximum  ice 
cover  The  work  was  done  as  part  of  a study  of 
nearshore  sea  ice  conducted  jointly  by  the 
United  States  Geological  Survey  (USGS)  and 
CRREL 

The  SLAR  used  was  a Motorola  APS-94-XF1 
which  operates  at  a frequency  of  10  GHz  The 
radar  was  mounted  below  a Mohawk  aircraft 
owned  by  the  USGS  The  ground  coverage  at  an 
altitude  of  1830  m was  in  two  swaths.  25  km 
wide,  on  both  sides  of  the  aircraft  A blind  spot 
10  km  wide  was  located  beneath  the  aircraft  and 
between  the  swaths 

The  SLAR  image  of  part  of  the  Coastal  Plain  is 
shown  in  Figure  23  An  examination  of  the  lake 
regions  reveals  that  some  lakes  give  strong 
returns  (bright  images)  while  others  give  weak 
returns  (dark  images)  Some  of  the  lakes  produce 
strong  returns  in  the  central  areas  but  give  weak 
returns  around  their  margins  An  interpretation 
of  this  difference  in  returns  (see  Weeks  et  al 
1977),  which  has  been  verified  by  limited  ground 
reconnaissance,  is  as  follows  The  strong  returns 
are  produced  at  the  strong  dielectric  permittivi- 
ty contrast  between  ice  (x  = 3 2)  and  water  (*'  = 
b0  at  100  MHz)  which  exists  in  the  deeper  lakes 
that  have  not  completely  frozen  The  weak 
returns  indicate  areas  where  no  significant  sub- 
surface dielectric  contrast  exists,  such  as  be- 
tween ice  and  frozen  soil  (4<k’<3)  The 
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a.  The  incident  signal  refracts  into  the  ice,  is 
backscattered  from  air  bubbles  and  then 
strongly  reflects  from  the  ice/water  interface. 


fa.  The  ice/ frozen  soil  interface  does  not  pro- 
vide a strong  dielectric  permittivity  contrast, 
allowing  the  signal  to  refract  strongly  into  the 
frozen  soil. 


Figure  24.  Possible  wave  reflection  and  refraction  mechanisms  according  to  Weeks  (1978)' 
that  may  explain  differences  between  radar  images  of  various  lakes  on  the  Arctic  Coastal 
Plain. 


mechanism  by  which  the  outgoing  signal  is 
backscattered  is  subject  to  debate,  but  Weeks* 
believes  it  may  result  from  the  numerous 
elongated  and  vertically  oriented  air  bubbles 
trapped  within  the  ice.  This  mechanism  is  il- 
lustrated in  Figure  24. 

Although  the  SLAR  imagery  appears  to  con- 
tain much  information  on  the  distribution  of 
lakes  with  water  available  under  the  ice  cover, 
there  is  no  way  in  which  to  determine  the  quanti- 
ty of  unfrozen  water  beneath  a lake  Therefore, 
once  the  appropriate  lakes  have  been  de- 
termined from  the  SLAR  imagery,  ground  recon- 
naissance must  be  used  to  determine  which  of 
these  lakes  is  best  suited  for  water  supply. 
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APPENDIX  A.  SATELLITE  IMACERY  FOR 
SUBSURFACE  WATER  EXPLORATION 


In  this  section  we  discuss  the  use  of  satellite 
imagery  at  optical  wavelengths  for  subsurface 
arctic  water  exploration.  Subsurface  water  is  in- 
ferred from  surface  features  which  indicate 
whether  or  not  a lake  is  totally  frozen.  Although 
optical  imaging  is  technically  an  electromag- 
netic technique,  it  is  discussed  in  an  appendix 
because,  in  contrast  to  most  other  geophysical 
techniques,  it  is  passive  and  has  marginal  earth 
penetration.  It  is  also  subject  to  weather  condi- 
tions. 

This  discussion  is  taken  from  a larger  study  by 
Sellmann  et  al.  (1975)  of  ERTS-1  (Earth  Resources 
Technology  Satellite)  imagery.  The  images  were 
obtained  throughout  the  spring  thaw  of  1973  at 
MSS  (multi-spectral  scan)  band  5 (0.6-0. 7 pm 
wavelength)  and  band  7 (0.8-1 .1  pm)  and  are  of 
the  same  Arctic  Coastal  Plain  Region  shown  in 
Figure  23.  The  features  of  interest  in  this  figure 
are  the  lakes  that  are  deep  enough  to  retain  sub- 
surface water  during  the  winter.  The  lakes  were 
divided  into  three  general  depth  categories 
based  on  changes  seen  in  sequential  imagery  ob- 
tained throughout  the  thaw  period  and  on  dif- 
ferences detected  between  bands  5 and  7.  The 
categories  are  1)  lakes  in  excess  of  2-m  depth 
that  do  not  freeze  to  the  bottom,  2)  shallow 
lakes  of  intermediate  depth  (1-2  m)  that  freeze  to 
the  bottom,  and  3)  very  shallow  lakes  and  ponds 
(<1  m),  some  of  which  are  basins  that  are 
seasonally  filled  at  the  time  of  spring  runoff. 


The  first  depth  class  (>2  m)  contains  the  im- 
portant lakes  in  terms  of  providing  winter  water 
supplies.  This  depth  class  can  be  detected  by 
observing  whether  or  not  a lake  in  question  has 
retained  its  ice  cover  long  after  the  shallow  lakes 
have  become  ice-free.  During  June  and  early 
July,  an  open  water  moat  corresponding  to  the 
shallow  water  areas  usually  forms  around  the  re- 
maining ice  of  the  deeper  central  basin  which 
can  be  as  much  as  2 m thick.  Figures  Al  (June) 
and  A2  (July)  at  band  7 show  that  the  decrease  in 
size  of  the  persisent  ice  masses  can  be  followed 
if  cloud  cover  permits  adequate  sequential 
coverage. 

The  second  depth  class  (1-2  m),  less  obvious 
than  the  first  group,  includes  shallow  lakes  that 
have  frozen  to  the  bottom,  usually  throughout 
their  basins.  In  Figure  A3  at  band  7,  these  lakes 
appear  dark,  suggesting  no  ice  cover  such  as 
noted  in  the  vicinity  of  Sinclair  Lake.  However, 
this  class  of  lakes  appears  white  at  band  5.  The 
difference  appears  to  be  caused  by  the  greater 
penetration  of  band  5 wavelengths  through  the 
ice  flooding  caused  by  surface  runoff  [penetra- 
tion of  light  through  water  for  band  5 radiation 
has  been  reported  to  be  as  great  as  3 m (Polcyn 
and  Lyzenga  1973)].  If  the  lakes  were  not  frozen 
to  the  bottom,  the  ice  mass  would  float  and  not 
be  subjected  to  surface  runoff,  thereby 
eliminating  the  contrast  betwen  bands  5 and  7. 

The  third  depth  class  (<1  m)  consists  of 
shallow  lakes  that  appear  dark  on  bands  5 and  7 
early  in  the  melt  season.  The  shallow  depths  of 
these  lakes  allow  the  ice  to  completely  melt 
quite  early  and  therefore  produce  no  contrasts 
between  the  different  wavelength  images. 
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A facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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